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I. INTRODUCTION 



When considering matter effects on neutrino oscillations, it is customary to consider 
only the charged current interaction of the electron-neutrino mediated by W^-exchange, and 
ignore the neutral current interactions of all three neutrino flavors mediated by Z-exchange. 
This is because the universality of the neutral current interaction ensures that the phases 
acquired by the three neutrino flavors through Z-exchange remain the same, and thereby 
do not lead to extra mixing effects beyond that due to W -exchange. 

However, in many models beyond the Standard Model (SM), the universality of the Zu e u t 
(£ = e, /i, r) couplings can be violated through radiative corrections, such as in Supersym- 
metric models with R-parity violating interactions [1], or through the mixing of the light 
active neutrinos with heavy sterile ones [2]. The existence of a Z' which couples to the 
three lepton flavors differently can also mimic the violation of universality in Z-exchange 
[3]. Though the violation of Zv{i>n coupling universality in the particular models considered 
in Refs. [1], [2], and [3] are strongly constrained by the universality of the Wlvi and Zii 
couplings, they nevertheless provide existence proofs that the universality of neutral current 
interactions cannot be taken for granted. 

The experimental bound on the violation of ZviVs coupling universality is also very weak. 
The sole constraint comes from CHARM and CHARM II [4, 5]: 

g Ue = 0.528 ± 0.085 , 
g v » = 0.502 ±0.017, 
g^/g"* = 1.051&H = 0.87-1.20, (1) 

where g ue is the coupling of neutrino flavor i/£ to the Z, normalized to 0.5 for the SM. These 
values were obtained from the measurements of the ratio and the double ratio R e / 
where 

_ ojyiN -> u t X) 
Rf = a(u e N =; l-X) ■ (2) 

The constraint on g"? was obtained from R^, and the constraint on the ratio g Ue / ' g Vti was 
obtained from the double ratio R e jR ll assuming charged current universality. The constraint 
on g Ve was obtained from those on g"* and g Ue /g Vfl . 

As we can see from the above numbers, while g v * is fairly well constrained to the SM 
value of 0.5, g Ve is ill constrained and can deviate significantly from 0.5. Of course, the sum 
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of squares of the ZvgVg couplings, namely 

{g^f + {g^f + GT) 2 , (3) 

is well constrained to its SM value by the Z invisible width measured by LEP and SLD 
[6], so any deviation in g" e must be accompanied by a corresponding deviation in g Vr to 
maintain this agreement. However, as long as g Ue and g Vr conspire to do so, large violations 
of universality are allowed. 

In this paper, we investigate the effect of such violations of neutral current universality 
on neutrino oscillations in matter. If the violation is as large as that allowed by CHARM 
and CHARM II, then it could lead to new effects that are measurable by long baseline 
neutrino oscillation experiments. If such effects are not seen, it could then improve upon 
the CHARM/CHARM II universality constraint. 

This paper is organized as follows: In section II, we derive the effective potentials due 
to the charged- and neutral-current interactions which enter the effective Hamiltonians that 
govern neutrino and anti-neutrino propagation in matter. In sections III and IV, we ap- 
proximately diagonalize the effective Hamiltonians for neutrino (III) and anti-neutrino (IV) 
propagation using the method of Ref. [7], and show how the effective mass-squared differ- 
ences and effective mixing angles are affected by the presence of neutral current universality 
violation. In particular, we will show that the effective mass-squared differences are little 
affected, while the shifts in the effective mixing angles are confined to just one angle; which 
angle this is depending on the mass hierarchy, and on whether the neutrino or anti-neutrino 
case is being considered. In section V, we discuss how these shifts in the effective mixing 
angles will manifest themselves in the neutrino and anti-neutrino oscillation probabilities, 
and point out that whether any effect can be seen or not depends crucially on the value 
of sin 2 (26*23) in vacuum. In section VI, we present the results of a numerical calculation of 
the effective mass-squared differences, effective mixing angles, and oscillation probabilities, 
which validate the approximations used in the previous sections. In section VII, we consider 
a hypothetical experiment in which the Fermilab NUMI beam [8, 9] in its high-energy mode 
is aimed at a 1 Megaton class detector 9120 km away at Kamioka, Japan (the planned 
Hyper-Kamiokande [10]) and discuss the potential constraint such an experiment can place 
on neutral current universality violation. Section VIII concludes. 
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II. THE EFFECTIVE POTENTIALS DUE TO W AND Z EXCHANGE 



Let us first derive the effective potentials for neutrino propagation in matter, which ac- 
count for the W- and Z-exchange interactions between the neutrino and the matter fermions. 
The effective potential due to ^-exchange is well known [11], but we will re-derive it in the 
following to provide a parallel to the Z-exchange case. 

At momentum transfers much lower than the W and Z masses, the weak interaction 
Hamiltonian of the neutrinos is given by 

-ffweak = H C c + H NC , (4) 
where Hcc and Hnc are the charged and neutral current contributions, respectively: 

h cc = 5| [p £ r (i - 75) e\ [H (1 - 75) M , 

G r - / \ 

Hnc = P -y| (! - 75) vt\ [/7 M [9v - 9a1s) f\ • ( 5 ) 

Here, £ is the lepton flavor {£ = e, /i, r); / denotes a generic fermion, and g v and are its 
vector, and axial-vector couplings to the Z: 

g f v = l(-2Q f sm 2 9 W: 

9 f A = 4 ■ (6) 

il and Qf are, respectively, the isospin and electric charge of the fermion /. If lepton-flavor- 
universality is violated in the neutral current interaction, the product pGp in the expression 
for if jvc wm depend on £. 

After a Fierz transformation, H C c can be rewritten as 

H cc = ^ W (1 - 75) €] [vav (1 - 75) ^] • (7) 

The forward scattering amplitude of a neutrino V£ against a non-relativistic lepton £ via 
W^-exchange is then 

Mcc = ^= (£\ £y° (1 - 7s) t \f) H mo (1 - 7s) v t h> 
= y/2 G F (£\ £H \£) (v t \ v\ (^^) v t \v t ) 



= V2G F N e (<f>l<j> ue ) , (8) 
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where N e = (£\ ftl \£) is the density of the charged lepton £, and <$> Vl is the two-component 
wave-function of the left-handed neutrino up. This shows that the effective potential that 
the neutrino experiences as it travels through matter is 

Vcc = V2G F N e . (9) 
In ordinary matter, iV M = N T = 0. Therefore, 

v cc M 2GfN ° (fOT " J ' do) 

y (for Vp, u T ) . 

Similarly, the forward scattering amplitude due to Z-exchange between a neutrino and a 
non-relativistic fermion / is given by 

G - / \ 
M NC = ~^j= (f\ /7° [9v ~ 9a1s) f \f) ( u e\ ^7o (1 - 7s) vt We) 

= \PlGp (f\ g f v (ff) \f) H v\ {^r) v t l^> 

= V2G F g f v N f {4>l e K) , (11) 

where we have set the p-parameter to one, and Nf = {f\ pf \ f) is the density of the fermion 
/. The effective potential due to the neutral current interaction is then 

V NC = V2G F g f v N f = V2G F (l( - 2Q f sin 2 6w) N f . (12) 

Since N e = N p in electrically neutral matter, we find 

Vnc — V2G F 

= - l - (V2G F N n ) . (13) 



\ + 2s 2 w \ N e +(\- 2^) N p+(-\) A '« 



Assuming N = N e = N p ^ N n , which is valid for the lighter nuclei which constitutes 
most of the Earth, we can relate N (cm' 3 ) to the matter density p (g/cm 3 ) via the Avogadro 
number N^: 

2N = N p + N n = pN A . (14) 

Then, 
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where we have used G F /(Hc) 3 = 1.16637(1) x 1(T 5 GeV~ 2 , N A = 6.0221415(10) x 10 23 mol -1 , 
and He = 0.197326968(17) GeV • fm [5]. Therefore, 

a . 2EV CC = (7.6324 x 10- eV 2 ) x (-^) x (^) , 

b = 2EV NC = ~a. (16) 

For anti-neutrinos, both Vcc and Vnc reverse their signs. 



III. THE EFFECTIVE MIXING ANGLES, NEUTRINO CASE 

A. Inclusion of Neutral Current Effects into the Effective Hamiltonian 

The effective potentials derived above enter the effective Hamiltonian for neutrino oscil- 
lations (multiplied by 2E) as follows: 
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Here, U is the MNS matrix in vacuum [12], a comes from the ^-exchange interaction of v e 
with the electrons in matter, while b e , 6 M , and b T come from the Z-exchange interaction of 
each neutrino flavor with the neutrons. If b e = b ll = b T = b, then the 6-matrix is proportional 
to the unit matrix, and it will not contribute to neutrino oscillations. However, if neutral 
current universality is broken, then b e ^ 6 M ^ b T in general and the 6-matrix cannot be 
ignored. 

The experimental constraints from CHARM/ CHARM II, Eq. (1), allow b e and b T to 
deviate significantly from b = —a/2, provided that b e + b T = 2b to satisfy the Z invisible 
width constraint, Eq. (3). We therefore write 

^ = l + 2£, ^ = 1, y = l-2£, (18) 
and use £ to parametrize the violation of universality. 
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Notice that we can rewrite the matter-effect terms in Eq. (17) in several different ways: 
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bu - b T 



+ 



6 M + b 7 



1 
1 
1 



• (19) 



Since the unit matrix terms can be dropped, this shows that we can always reduce the 
problem to the case b T = 0, or 6^ = —b T . We will use the latter replacement in the following. 
For the case of Eq. (18), this entails making the replacement 



a + 6(1 + 20 
6 

6(1 - 20 



(a + 360 
6£ 
-&£ 



(20) 
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a 





H = U 
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5m 2 31 










Furthermore, we absorb the factor 36^ in the (1, 1) element into a since we can expect 
36£ a, and the uncertainty in the matter density p which enters into a can be expected 
to hide any such shift. Therefore, the effective Hamiltonian we will consider is 



(21) 



The problem is to diagonalize H and find the eigenvalues A, (i = 1,2,3) and the diagonal- 
ization matrix U. 

To this end, we use the method of Ref. [7] in which the Aj's and U were derived for 
the £ = case. The procedure followed in Ref. [7] was to approximately diagonalize the 
effective Hamiltonian, H, using the Jacobi method: 2x2 submatrices of H are diagonalized 
in the order which requires the the largest rotation angles until the off-diagonal elements are 
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negligibly small. As the order parameter to evaluate the size of these off-diagonal elements, 
we use 



Sm 2 
1 5m 2 



21 



(22) 

<-3ii 

and consider H to be approximately diagonalized when the rotation angles required for 
further diagonalization are of order e 3 or smaller. For Sm^ = 8. 21°;!? x 10" 5 eV 2 and \5m 31 | = 
(1.5 ~ 3.4) x 10" 3 eV 2 , we have e = 0.15 ~ 0.24 and e 3 = 0.0034 ~ 0.014. 

For the sizes of the mixing angles in vacuum, we assume 9 13 = 0(e), cos(2# 12 )/2 = 0(e), 
and cos(2# 23 ) < 0(e) as in Ref. [7]. We also assume that the universality violation parameter 
f is of order e 2 = 0.02 ~ 0.06, since the central value of the CHARM/CHARM II result 
translates to £ = 0.025. 



B. Diagonalization of the Effective Hamiltonian 



For the neutral current term in Eq. (21) to have a non- negligible effect on neutrino 
oscillations, we anticipate that it must be at least as large as, or larger than, the smaller 
mass-squared-difference Sm^. Since we have assumed £ = 0(e 2 ), this requires a = —2b = 
2\/2G F NE to be at least as large as, or larger than, the larger mass-squared-difference 
Itfmfj-J. For the sake of concreteness, we will consider the case a/\8m\ 1 \ = 0(e~ l ) in the 
following. 

Introducing the matrix 

Q = diag(l,l,e w ), (23) 



we begin by partially diagonalizing the Hamiltonian H as 



H' = Q)U ] HUQ 





5m 2 21 
5m 2 31 



5m 2 21 
Sm 2 31 



a 
K 
-&£ 



UQ 



u* el u el u* el u e2 u* el u e3 
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U* e3 U el U* e3 U e2 U* e3 U e3 
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u;^ u;,u, 2 u; x u^ 
+bi Qt { u* 2Ufxi 

u; 3 u^ u; 3 u, 2 u; 3 u, 3 

The matrix multiplying a is given by 



U* Tl U Tl U* t1 U t2 U* t1 U t3 

u; 2 u Tl u; 2 u r2 u; 2 u t3 
u; 3 u Tl u* t3 u t2 u; 3 u t3 



)Q. (24) 



M a = Q 



u* el u el u* el u e2 u* el u e3 

U* e2 U el U* e2 U e2 U* e2 U e3 
U* e3 U el U* e3 U e2 U* e3 U e3 



Q = 



r 2 r 2 

c 12 c 13 



Cl2<Sl2 c i3 



Cl 2 S\ 2 C l3 C12C13S13 



„2 2 
*12 c 13 



C12C13S13 S12C13S13 



512C13S13 



'13 



Using 6*13 = 0(e), we estimate the sizes of the elements of M a to be 



M„ 



0(1) 0(1) 0(e) 
0(1) 0(1) 0(e) 
0(e) 0(e) 0(e 2 ) 



The matrix multiplying is given by 

Uk u i* U U U ^ U U U ^ 



M b = QU 



u; 2 u»r u; 2 u, 2 u; 2 u, 3 
u; 3 u^ u; 3 u, 2 u; 3 u, 3 



u* Tl u Tl u* Tl u T2 u* t1 u t3 
u: 2 u Tl u; 2 u r2 u* t2 u t3 
u: 3 u Tl u; 3 u r2 u; 3 u t3 



> Q 



(s\ 2 - c\ 2 s\ 3 ) cos(20 23 ) + sm(20 12 ) sm(2d 23 )s 13 cos5 
-(1 + s 2 3 )s 12 c 12 cos(20 23 ) - (cj 2 e- lS - s 2 2 e +i5 )s 13 sm(29 23 ) 
-S12C13 sin(29 23 )e~ lS + ci 2 s 13 c 13 cos(29 23 ) 

-(1 + s 2 3 )s 12 c 12 cos(29 23 ) - (c 2 2 e +lS - s 2 2 e- i5 )s 13 sin(20 23 ) 



-12 



„2 2 



> 12 ° 13 ) cos(26 23 ) - sin(2#i 2 ) sin(2^23)si3 cos S 
c 12 c i3 sm(29 23 )e~ lS + s 12 s 13 c 13 cos(26> 23 ) 

-S12C13 sin ( 26*23 )e +i5 + c 12 s 13 c 13 cos(29 23 ] 
c 12 c 13 sin ( 26*23 )e +i<5 + S12S13C13 cos(26 23 ] 
-c\ 3 cos(29 23 ) 



(25) 



(26) 



(27) 



Using cos(29 23 ) < 0(e) and 6*13 = 0(e), we estimate the sizes of the elements of M b to be 



0(e) 0(e) 0(1) 
0(e) 0(e) 0(1) 
0(1) 0(1) 0(e) 



(28) 
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Since we are only interested in the leading order effect in £, we neglect the 0(e) terms in 
Mb and approximate 













-s 12 e tS c l2 e tS 



si 2 e 
cue 




Under this approximation, the effective Hamiltonian that must be diagonalized is 



H' = diag(0, 5m 2 21) 5m 2 31 ) + aM a + b£M b 



2 2 2 
aC 12 C 13 fl Ci2Si2C 13 

a C12S12C13 a s\ 2 c\ 3 + 5m 2 .! 

a c 12 c l3 s 13 - b£ s 12 e~ tS a s l2 c l3 s 13 + b£ c X2 e 



ac 12 c 13 s 13 - b£s 12 e +tS 
as 12 c 13 s 13 + b£ c 12 e +tS 



a s; 



+ 5ml 



At this point, we set 5 = for the sake of simplicity. Then H' becomes 



H' = 



a C 12 C 13 



ac 12 si 2 c 13 



a c 12 s 12 c 2 13 a c u c 13 s 13 - b£ s 12 
a s\ 2 c\ 3 + 5m 2 21 a s 12 ci 3 s 13 + b£ c 12 

is + 5m ti 



aci 2 ci 3 si 3 - b£s 12 asi 2 c 13 sis + b^ci 2 as 2 ^ + 5m 2 



(29) 



(30) 



0(1) 0(1) 0(e) 
0(1) 0(1) 0(e) 
0(e) 0(e) 0(e) 
(31) 



Bl. First Rotation 



Applying the Jacobi method to H', we first diagonalize the (1, 2) submatrix which requires 
a rotation by an angle of 0(1). Define the matrix V as: 



Then, 



H" = V ] H'V = 












V = 


Sip C(p 





(32) 







1 





where 



dn = cos ip , s v = sin ip , tan 2ip 



ac 



\ 3 sin 26 



12 



5m 2l — ac\ 3 cos 26> 12 



(O < <p < \ ~ On) ■ (33) 











ac' 12 C13S13 - fc£s' 12 
asi 2 ci 3 si3 + &£c' 12 



ac' 12 ci 3 si3 - 6£s' 12 asi 2 ci 3 si 3 + b£c' l2 as\ 3 + 



(34) 
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where 



and 



c'i2 = cos 6' 12 , s' 12 = sin 6[ 2 , 9' 12 = 9 12 + (p , 



(ac? 2 c2 3 ) C 2 - (as 2 2 c 2 3 + ^mljsj 
A i ~~ 3 ~ ~~ A - > 



c% - si 



a; 



(«5? 2 c? 3 + bm 2 21 )c% - (ac 2 2 c 2 3 )s 2 

- A + > 



c 1 s 1 



with 



a: 



(ac 2 3 + o~m 21 ) ± yj (acf 3 — 5m 2l ) 2 + 4ac 2 3 sf 2 o~m 



21 



^-^isin(2M + 0(. 5 ), 



As discussed in Ref. [7], in the region a/lom^l = 0(e 1 ), we can expand 9' 12 as 
from which we can conclude 

1.?. 

■ sin(2# 12 



S 12 
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Also, expanding A' ± , we find 
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find 



cos 



sin(2tf 



sin(26 l i2) I = sin 



2a 
2a 



12, 



sin(2# 12 ) 



1 - 0(e & ) 
0(e 3 ) . 



A'_ = 5m^4 + aO(e 6 ) = aO(e 3 ) , 

A' + = ac 2 3 + ^ m 2i s i2 + aO(e 6 ) = a 0(1) . 



Therefore, the sizes of the elements of H" are evaluated to be 



H" 



0(e 3 ) 0{e 2 ) 
0(1) 0(e) 
0(e 2 ) 0(e) 0(e) 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



Unlike the £ = case considered in Ref. [7], both the (1,3) and (2,3) submatrices require 
rotations by angles of 0(e) to diagonalize. Here, we diagonalize the (2,3) submatrix next 
to maintain the parallel with the £ = case. 
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B2. Second Rotation 



The matrix W which diagonalizes the (2, 3) submatrix is 



W 



1 


-s^ 



(42) 



where c<f, = cos0, s<j, = sin0, and 



, o , 2(as' 12 s 13 c n + b£c' 12 ) 

tan 20 = — - — ^ ^ — — 

5m^ + as 2 13 - A' 2 



2a£i2gi3Ci3 

5ml + as?, — An 



"12 



2s' 12 Si3Ci3 



tan20 o <!l-q^T^ ) 



(43) 



The angle is in the first quadrant when 5m 31 > (normal hierarchy), and in the fourth 
quadrant when 5m 3l < (inverted hierarchy). O is the rotation angle when £ = 0. Taking 
the arc-tangent of both sides of the above equation, we find 



When a I \ 5m 



0(e 



4 sin 2#i3 
b is given by [7] 



sin 40o cot 0' 12 \ +0 (g2j 



(44) 



'13 



5m 



31 



— (713 



a 

0i 3 + 0(r 



^13 + 0(£ 3 ) 
3^ 



W > o) , 
W < o) • 



(45) 



Therefore, sin40 o 
Eq. (38), we find 



sin(40 13 ) for both the 5m^ > and (^m^ < cases and using 



~ 00 + £ 



^21 

4a 



sin(2# 



12 



(46) 



Therefore, the difference between and O can be neglected in this range. 
Using W, we obtain 



H'" = W^H"W 



A'i 



-(aci 2 ci 3 si 3 - 6f s' 12 )s0 (ac' 12 ci 3 si3 - 6£si 2 )c 



-(aci 2 ci 3 si 3 - 6£s 

12)^ 

(ac' 12 ci 3 s 13 - &£s' 12 )c 



A" 



A 3 



(47) 
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where 



A 2 — 



ci - si 



(as 2 13 + 5m 2 31 )c 2 - X' 2 s 2 



K = 



& - si 



(48) 



If we define 

.// _ t A 2 + (as? 3 + $ m li)} ± V^2 - (« s ?3 + Sm 2 31 )} 2 + 4(as / 12 c 13 s 13 + fefci 



A 



12 J 



± 



(49) 



then 



as = a': 



A 3 ~ A + > 



A 2 



A + ' A 3 — A - ' 

When a/l&njjjj = 0(e _1 ), we can expand A± as 



if 5m\ r > , 
if Sm^ < . 



(50) 



A" = 5m 2 31 c\ 3 + 0(e 3 \5m 2 31 \) = a 0(e) , 

A'| = a + 5m 2 31 s 2 13 + 5m 2 21 s 2 12 + 0(e 3 \5m 2 31 \) = a 0(1) 

Also, from Eq. (45) and the fact that (j) o , we conclude 



(51) 



si 3 = 0(e) , (5m 2 ! > 0) , 
cia = 0(1) , (8m 2 ! < 0) . 



(52) 



Sep ~ C13 = 0(1) , C0 

s<f> « -S13 = 0(e) , 

In the £ = case considered in Ref . [7] , H'" was already approximately diagonal and further 
diagonalization was not necessary. However, when £ = 0(e 2 ), the sizes of the elements of 
H'" are found to be 



H" 



0(e 3 ) 0(e 2 ) 0(e 3 ) 



0(e 2 ) 0(e) 
0(e 3 ) 





0(1) 



(53) 



when Sml ± > 0, and 



H" 



0(e 3 ) 0(e 3 ) 0(e 2 





0(e) 



(54) 



0(e 3 ) 0(1) 
0(e 2 ) 

when Sm^ < 0. Therefore, for the 5m 2 31 > case (normal hierarchy) we must diagonalize the 
(1, 2) submatrix next, while for the 5m 31 < case (inverted hierarchy) we must diagonalize 
the (1,3) submatrix next. 
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B3. Third Rotation, Sm^ > Case 



Define the matrix X as 



X 



where 



Then, 



c x s x 

~ s x c x 
1 



c x = cos x , s x = sin x , tan 2x 



-2(ac' 12 c l3 s 13 - 6£s' 12 )s0 



TTllll 



= X^H'"X 



A ix 



A 2X 



(ac' 12 c 13 s 13 - b^s' 12 )c^c x 
{ac' l2 c 13 si 3 - b£s' 12 )c <j) s x 



(ac' 12 c 13 s 13 — 6^s' 12 )c0C x (ac' 12 Ci 3 si3 — b^s'^c^s^ 



A 3 



where 



A ix — 



A 



2X 



A i c x 


X" e 2 
~ A 2 S X 


c 2 
°x 

\"r 2 


- S 2 

X 

\> « 2 


C 2 


- s 2 

b x 



— A x- > 



X'" 

A X+ > 



with 



A x± 



(K + A'x) ± \J (A 2 - Ai) 2 + 4(ac' 12 c 13Sl 3 - &f s' 12 ) V 



Recalling that 



ac 'l2 C 13 s 13 


/ 5ml, 
V 2a 


sin(26»i 2 )^ C13S13 = aO(e 4 ) , 




a 

2^ 


-aO(e 2 ) , 


a; 


~ 5m 21 c 12 


= aO{e 3 ) , 


A' 2 ' 


~ ^ m 31 c 13 


= aO(e) , 



and s«£ ~ 1, we find 



tan2 X = - 2 «y-^)^ = {1 + 0(. 2 )} 



A 2 - X[ 



5m 2 31 
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Therefore, the angle x is given approximately by 



X 



5m 2 31 



0(e), 



(62) 



from which we can conclude that s x = 0(e) and c x = 0(1). The eigenvalues can also be 
expanded in e and we find 



A'" 



IX 



A'" 



2X 



5m 2 21 c 2 12 - -\- + aO(e 5 ) = aO(e 3 ) , 
om 31 

b 2 £ 2 



+ aO(e A ) = a 0(e) . 



(63) 



Note that these shifts of the eigenvalues are of order aO(e 3 ) and have negligible effect on 
5\si = aO(l) or 6X21 = aO(e). Putting everything together, we evaluate the sizes of the 
elements of H™' to find 

r 0(e 3 ) 0(e 3 ) 

0(e) 0(e 4 ) ■ (64) 
0(e 3 ) 0(e A ) 0(1) 

This shows that further diagonalization requires rotations by angles of 0(e 3 ) or smaller, 
which we will neglect. 

Thus, we have found that when £ = 0(e 2 ) and Sm 31 > (normal hierarchy), we need an 
extra (1, 2)-rotation to diagonalize H, and the diagonalization matrix is UVWX, which we 
need to identify with 



U = 



1 

C23 S23 

— S23 c 23 



Cl3 

1 

id 



s^e~ iS 



C12 S12 
-S12 C12 
1 



(65) 





-si 3 e' u U ci 3 

to obtain the effective mixing angles and effective CP phase. From Ref. [7], we know that 
when fonij;! > anc ^ a /l^ m iil — 0( £-1 )j identification of U with UVW leads to 

$13 — 
012 



'13 ) 



2 

#23 — #23 

8 = 0, 

where we have defined 

^i 3 = ^13 + 0, 







L 13 v 


v 2a y 






c 13 


V 2a 



sin(2fl 12 ) , 
sin(2# 12 ) 



-13 



COS 0' , 



s' 13 = sin ^ 3 . 



(66) 



(67) 
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(Recall that we are considering the 5 = case). Since X is an (1, 2)-rotation, multiplication 
of U — UVW from the right by X only shifts the value of 6>i 2 by \- Therefore, we can 
conclude that 



#13 
#12 

#23 

6 



13 ) 
7T C13 



-13 



#23 + 







13 



^21 

2a 
2a 



sin(2# 12 ) + x , 
sin(2fl 12 ) , 



(68) 



In these expressions, non only x, but also <p and 9' 13 = 9 13 + <f> depend on £. However, the 
^-dependence of is very weak. The ^-dependence of 5\ 3i and 5X 2 i are also weak, so the 
effect of a non-zero £ will appear dominantly in #12. 

.EL£. Third Rotation, 5m 31 < Case 

In this case, we need to diagonalize the (1, 3)-submatrix of H'". Define the matrix Y as 



Y = 



1 



Sfj Cq 



where 



Then, 



c v = cost] , s v = sin?7 , tan2r? 



2(ac' 12 c 13 s 13 - b^s'^c^ 
A 3 - X 1 



(69) 



(70) 



Ky -(ac' 12 ci 3 s 13 - 6£s' 12 )s c^ 

-(aci 2 ci 3 si3 - 6£s' 12 )s0C^ A 2 

-(aci 2 ci 3 si 3 - b^s'^s^ 



[ac' l2 ci 3 s 13 - 6£s' 12 )s0S ri 



V" 



(71) 



where, 



x/// A lS _ x/// 

Aiv ~ ,2 _ e 2 ~~ A Y+ ' 



- S< 
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\"r 2 — \' q 2 
w// _ A 3 C rj A l b r t _ */// 

A 3Y - .2 _ e 2 ~ A Y- 



c v - 



with 



Using 



A Y± 



(AS + AD ± - Ai) 2 + 4(ac' 12 c 13S i3 + &^)^ 



(72) 



(73) 



ac' 12 c 13 s 13 


( 5ml, 

ks a — 

V 2a 




a 

i~*~> i- 

2^ 


K 


~ ^ m 21 C 12 


K 


~ ^ m 31 C 13 



sm(26»i 2 )J C13S13 = aO(e 4 ) , 

-aO(e 2 ) , 
= aO{e 3 ) , 
= a 0(e) , 



(74) 



and Ctf) ~ 1, we find 



tan2, = 2 «^i3 ~ *^ C ' = {1 + 0(e 2 )} 

A3 — A^ dm 



Therefore, the angle 77 is given approximately by 

bi bi 



V 



0(e), 



(75) 



(76) 



from which we can conclude that s v = 0(e) and c v = 0(1). The eigenvalues can also be 
expanded in e and we find 



Ky = <^21 C 12 



5m 2 31 



+ aO(e b ) = aO(e A ) , 



\& = 5mlA, + ^r + aO(e A ) = a 0(e) 



(77) 



Again, these shifts in the eigenvalues are negligible. Putting everything together, we evaluate 



the sizes of the elements of H'y to find 



Hy = a 



0(e 3 ) 0(e 3 ) 
0(e 3 ) 0(1) 0(e A ) 
0(e 4 ) 0(e) 



(78) 



This shows that further diagonalization requires rotations by angles of 0(e 3 ) or smaller, 
which we will neglect. 
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Thus, we have found that when £ = 0(e 2 ) and Sm^ < (inverted hierarchy), we need an 
extra (1, 3)-rotation to diagonalize H, and the diagonalization matrix is UVWY. We need 
to identify this product with 



U = 



1 

c 23 S 2 3 
-S 2 3 C 23 



ci 3 s i3 e 
1 



-j<5 



Cl2 S12 
-Sl2 C12 
1 



(79) 



-si 3 e*" U ci 3 

to obtain the effective mixing angles and effective CP phase. From Ref. [7], we know that 
when Smli < and a/l&njjjj = 0(£ _1 ), identification of U with UVW leads to 



Furthermore, in the range a/l&nij-J = 0(e x ) we have 







13 — 




12 — 




23 = 


#23 


5 = 


. 



H sin(2M + - = J + 0(£ 3 ), 



2a 



#13 - 



31 



#13 + 



which implies that 





1 










1 










1 





UVW « 





c 23 


s 23 







1 







-1 










~S23 


C 2 3 










1 







1 



Then, 



UVWY 



1 

c 23 S 2 3 

-S 2 3 C 23 

1 

U C 23 S23 

-s 2 3 c 23 



1 

1 

1 

1 



Cjj S-q 

Sjj 



1 
-10 
1 

1 
1 
1 



Cjj S-q 



1 



s r i 



1 
-10 
1 



Therefore, 77 can be absorbed into 6> 2 3 as 

013 = #13 > 
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(80) 



(81) 



(82) 



(83) 



012 — #12 > 
023 = 023 — V , 

8 = 0. (84) 

As in the 8rn^ 1 > case, the ^-dependence of 9' 13 = 0i3 + is very weak, so the effect of a 
non-zero £ will appear dominantly in #23- 



C. Summary of Neutrino Case 



To summarize what we have learned, the main effect of including the terms, which 
come from neutral current universality violation, in the effective Hamiltonian is to shift 9 12 
in the Sm^ > (normal hierarchy) case, and #23 in the Sm^ < (inverted hierarchy) case, 
beyond the shifts due to the charged current interaction term a. In the Sm 31 > case, the 
shift in 9 12 is given by 

X ~ = , (85) 



5m 2 31 



2 Smli 



while for the Sm 31 < case, the shift in #23 is given by 



-77 



^ m ii 



6£ 



2 1 I 



(86) 



IV. THE EFFECTIVE MIXING ANGLES, ANTI-NEUTRINO CASE 



Inclusion of Neutral Current Effects into the Effective Hamiltonian 



For the anti-neutrinos, the effective Hamiltonian is given by 





Ai 





















—a 










-be 








H = U* 





A 2 





U T = u* 





5m 2 21 





U T + 











+ 





-b. 













A3 










5m 2 31 
























(87) 

The differences from the neutrino case are the reversal of signs of the CP violating phase 5 
(and thus the complex conjugation of the MNS matrix U), and the matter interaction terms 
a, and bg (£ = e, //, r). We denote the matter effect corrected diagonalization matrix as U 
(note the mirror image tilde on top) to distinguish it from that for the neutrinos. As in the 
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neutrino case, we make the replacement 



—a 










-be 

















+ 





-K 






















-b T 



-a 
-&£ 
b£ 



(88) 



B. Diagonalization of the Effective Hamiltonian 

Using the matrix Q from Eq. (23), we begin by partially diagonalize the effective Hamil- 
tonian as 



H' 



QU T HU* Q* 




5777,2! 








5m 2 
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5m 2 31 




5m 2 31 



- QU 1 



aQ 



UplU*! U^U* 2 U^U* 3 

u, 2 u;, u, 2 u; 2 u, 2 u; 3 



UvaUpi UijzU* 2 U^U^ 3 



a 

hi V Q 
-&£ 

U el U* el U el U* e2 U el U* e3 
U e2 U* el U e2 U* e2 U e2 U* e3 Q 
U e3 U* el U e3 U* e2 u e3 u* e3 

u Tl u* Tl u Tl u; 2 u Tl u* r3 
u T2 u* Tl u t2 u; 2 u t2 u* t3 
u t3 u* t1 u t3 u; 2 u t3 u; 3 



> Q* . (89) 



The matrix which multiplies a is given by 

u el u* el u el u: 2 u el u* e3 
•u, = Q u e2 u* el u e2 u* e2 u e2 u* e3 q: 

U e3 U* el Ue3U* e2 u e3 u; 3 
while the matrix which multiplies is given by 



2 2 2 
C 12 C 13 c 12S\ 2 C l3 C12C13S13 

C12S12C13 s\ 2 C 2 l3 S12C13S13 



UniU*! 17^17*2 C/ M i[/* 3 

u, 2 u; x u, 2 u; 2 u, 2 u; 3 
u, 3 u;, u, 3 u; 2 u, 3 u; 3 



C12CI13S13 S12C13S13 



u Tl u;, u Tl u: 2 u Tl u; 3 

u r2 u* Tl U t2 W t2 U t2 U* t3 

u t3 u* t1 u t3 u; 2 u t3 u* t3 



'13 



} Q* 



(90) 



Mi 



( s i2 _ c i2 s i3) cos(2^23) + sin(2#i 2 ) sin(2^23)si3 cos 5 
-(1 + s 2 13 )s 12 c 12 cos(29 23 ) - (c 2 12 e +iS - s 2 2 e~ i5 )s 13 sin(20 23 ) 
-s 12 c 13 sin ( 26*23 )e +i<5 + c 12 s 13 c 13 cos(26> 23 ) 
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-(1 + s 2 13 )s l2 c 12 cos(2# 23 ) - (c\ 2 e- lS - s 2 2 e +t8 )s 13 sm(26 23 ) 
( c i2 - s i2 s i3) cos(26»23) - sm(29 12 ) sm(29 23 )s 13 cos5 
c 12 c 13 sin ( 26> 23 )e +lS + S12S13C13 cos(20 23 ) 

-S12C13 sin(26»23)e~ i ' 5 + c 12 s 13 c 13 cos(26> 23 ) 
C12C13 sin(26> 2 3)e^' 5 + s 12 s 13 c 13 cos(26> 23 ) 
-c\ 3 cos(2# 23 ) 



(91) 



Using cos(26> 23 ) < 0(e) and 6> 13 = 0(e), the sizes of the elements of M b are evaluated to be 



M h 



0(e) 0(e) 0(1) 
0(e) 0(e) 0(1) 
0(1) 0(1) 0(e) 



As in the neutrino case, we neglect the 0(e) terms in Mh and approximate 



M h 








-s 12 e 



-s l2 e +iS c 12 e +iS 



ci 2 e 




-iS 



The effective Hamiltonian which must be diagonalized is then 



(92) 



(93) 



H' = diag(0, 5m 2 21 , 5m 2 31 ) - a M a - b£ M b 

„2 „2 



a C 12 C 13 



-aci 2 si 2 c 13 
-a Ci 2 Si 2 cf 3 -a s\ 2 c\ 3 + tfm^ 

-a C12S13S13 + K s X2 e +lS -a s 12 c n s 13 - b£ c 12 e +tS 



-aci 2 Ci 3 s 13 + b£ si 2 e 
-as 12 c 13 s 13 - b£,ci 2 e~ 



—a s 



13 



5m 2 31 



From this point on, we set 5 = for the sake of simplicity. Then, H' becomes 



H' = 



-a c 2 c 2 



12^13 



-a c 12 s 12 cj 3 -a C12C13S13 + b£ s 12 
-a c 12 s 12 cj 3 -a s\ 2 c\ 3 + 5m 2 21 -a s 12 c 13 s 13 - b£ c 12 

-ac 12 s 13 s 13 + b£s 12 -a s 12 c 13 s 13 - b£ c 12 -as 2 13 + Sm 2 31 



= a 



(94) 



0(1) 0(1) 0(e) 

0(1) 0(1) 0(e) 

0(e) 0(e) 0(e) 
(95) 
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Bl. First Rotation 



Applying the Jacobi method on H', we begin by diagonalizing the (1, 2)-submatrix. Define 
the matrix V as 



V = 

where 

c v = cos <f , s v = sin <p , tan 2<p = 
Then, 



1 



ac 2 3 sin 2# 12 



5m 2l + ac 2 3 cos 29 



12 



if" = 0ff' y 







An 



-ac' 12 c 13 s 13 + 6£s' 12 
-as' 12 c 13 s 13 - b£c' 12 



where 



and 



-ac' 12 c 13 s 13 + b£s' 12 -as' 12 c 13 s 13 - b£c' 12 -asj 3 + 5m 2 31 



C'i2 = COS Q'l2 » S 'l2 = Sm 0l2 > 'l2 = ^12 + V , 



Ai = 



K = 



-2 „2 



{-ac{ 2 c{ 3 )c' + {as\ 2 c{ 3 



$rn 2 2l )sl 

— A - , 



as\ 2 c\ 3 + om^c 2 , + (ac 2 9 c 2 3 )sl 



°21 J ^ifi 



12 < -'13,' C V 



with 



A' 



(5m 21 — ac 2 3 ) ± a/ (5m 2l + ac 2 3 ) 2 — Aac\ 3 s\ 2 5m 



21 



± 



From Ref. [7], we know that in the region a/l&njj-J = 0(e 1 ), we can expand (p as 



^2 = ^r i sin(2M + 0(e 5 ) = 0( £ 3 ) 
2a 



Therefore, 



C 12 



S 12 



COS 



sm 



^21 

2a 
<5m 2 



sin(2£ 12 ) =l-0(e 6 ), 



21 



2a 



sin(2fl 12 ) = 0(£ 3 ) 



(96) 



(-012 < V? < 0) . (97) 



(98) 



(99) 



(100) 



(101) 



(102) 



(103) 
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The expansions of A' ± are given by 



A' 



— ac^ 3 + 5m 2 A s 



2i°i2 + a O(e t 



a 0(1) , 



\' + = 5m 2 21 c 2 12 + aO(e 6 ) = aO(e 3 ) . 
Therefore, the sizes of the elements of H" can be evaluated to be 

0(1) 0(e) 

H" 



(104) 



(105) 



0(e 3 ) 0(e 2 ) 
0(e) 0(e 2 ) 0(e) 

As in the neutrino case, though we have a choice of whether we diagonalize the (1,3) or the 
(2,3) submatrix, since both require rotations by angles of 0(e), we diagonalize the (1,3) 
submatrix next to maintain the parallel with the £ = case. 



B2. Second Rotation 



Define the matrix W as 



W = 



c<p s<f> 
1 
-s^ 



(106) 



where = cos0, = sin0, and 



tan 20 = -2(^13** - 6^) 



5m 31 — as 2 3 — X[ 
2ac 12 Ci3Si3 



Sm^ — as\ z — X[ 



1 + C 



'12 



2c' 12 Ci3Si3 



= tan I 



1+e 



'12 



(107) 



2c' 12 Ci3Si3 / 

The angle is in the fourth quadrant when 5m 31 > (normal hierarchy), and in the first 
quadrant when 5m 31 < (inverted hierarchy). 0o is the rotation angle when £ = 0. Taking 
the arc-tangent of both sides of the above equation, we find 

'sin 40o tan 9' 12 " 



= 00 + £ 



When a/l^mU = 0(e O is given by 



4 sin 29 



13 



+ o(r 



(108) 



00 



'13 



Sm 2 31 



'13 

<5m 



(K > o) 



(109) 
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Therefore, sin40 o ~ — sin(4#i 3 ) for both the Sm^ > and Sm^ < cases and using 
Eq. (102), we find 

4> « 0o - $ ( ^fO sin(2# 12 ) = O - ^(e 5 ) . (110) 

Therefore, the difference between and 0o can be neglected in this range of a, just as in the 
neutrino case. 

Using W, we obtain 



H'" = W ] H"W 



(as' 12 s 13 c 13 + b£d 12 )s4 



(as' 12 s 13 c 13 + &£c' 12 )s 

K 

-(as' 12 s 13 c 13 + btc'^cj, 





-(as' 12 s 13 c 13 + ftfcj 



12) 



(AH) 



where 



K = 



K4 + (ash ~ Sm 2 31 )sl 



r 2 - s 2 



(-as 2 3 + (^mlJc 2 - A' lS 2 



r 2 - s 2 



(112) 



If we define 



// _ i( Sm ii - as is) + K) ± V#4 - as ?s) - A iP + 4:(ac[ 2 s 13 c 13 + b£s' 12 ) 2 



(113) 



then 



A'/ = A" , A 3 ' = A'j , if 5m 2 3l > , 
A" = A'| , A3 = A" , if 5m 2 1 < . 



(114) 



When a/l^mU = 0(e 1 ), we can expand A± as 



A" = -a + s 2 13 Sm 2 31 + s\ 2 8m 2 21 + 0(e 3 \Sm 2 31 \) = aO(l) , 



A'| = c? 3 5m^ + 0(e 3 |<5m 



311 



a 0(e) . 



(115) 



Also, from Eq. (109) and the fact that 



we conclude 



S4, ~ -S13 = O(e) , « C13 = 0(1) , (5m| x > 0) , 
~ ci 3 = 0(1) , c,,-, « si 3 = O(e) , (5m| x < 0) . 



(116) 
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Putting everything together, we evaluate the sizes of the elements of H'" and find 



when Sm 31 > 0, and 



H'" = a 



0(1) 0(e 3 ) 
0(e 3 ) 0(e 3 ) 0(e 2 ) 
0(e 2 ) 0(e) 



(117) 



0(e) 0(e 2 ) 
//"' a 0(e 2 ) Oi:-) 0(e 3 ) (118) 
0(e 3 ) 0(1) 

when Smli < 0. Therefore, for the Sm^ > (normal hierarchy) case, we must diagonal- 
ize the (2, 3) sub matrix next, while for the 5m 31 < (inverted hierarchy) case, we must 
diagonalize the (1,2) sub matrix next. 



B3. Third Rotation, 5m 31 > Case 



To diagonalize the (2, 3) submatrix of H'" , we define the matrix X as 



X 



1 



c x s x 
~s x c x 



where 



Then, 



c x = cos x , s x = sin x , tan 2% 



A 3 ~~ A 2 



TTllll 

H x = 



(119) 



(120) 



\'( (as' 12 s 13 c 13 + 6^c' 12 )s^c x (as' 12 s 13 c 13 + b^c' 12 )s 4> s x 

as' 12 s 13 c 13 + ^ci 2 )^c x \' 2 ' x 

as' 12 s 13 c 13 + b£c' l2 )s ( j > s x \% x 



{121) 



where 



Xlll _ A 2 C X - A 3 S X 
A 2X - 



C X S X 



A'" 



X- i 
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\//-2 \l ^2 
A 3X — — A 



x+ > 



with 



Recalling that 



X w — 
A x± = 



(AS + A' 2 ) ± - A' 2 )2 + 4(as' 12Sl3 c 13 + &e^ 2 ) 2 ^ 



122) 



(123) 



as' 12 s 13 c 13 p 


S ° { 2a 




2^ 




s 5m 21 c 12 


A 3 " 





sin(20i 2 ) ) S13C13 = «0(e ) , 
-aO(e 2 ) 



(124) 



and Ctj, ~ 1, we find 



2 fas 7 . 



tan2x = - -^-^13013 + ^)^ = _ 26|_ {l Q 2 } 



A 3 ~~ A' 2 



(125) 
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Therefore, the angle x is given approximately by 



X 



= 0(e) , 



(126) 



from which we can conclude that s x = 0(e) and c x = 0(1). The eigenvalues can also be 
expanded in e and we find 

b 2 e 



X'2'x = Sml^c 2 



21^12 



5m 



+ aO(e 5 ) = aO(e 3 ) , 



31 
2f2 



K'x = W^ + ^ + aO^ 4 ) = a 0(e). 



(127) 
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As in the neutrino case, the shifts in the eigenvalues are of order aO(e 3 ) and their effects 
on <5A 3 i = aO(l) and 5A 2 i = aO(l) are negligible. Putting everything together, we evaluate 



the sizes of the elements of H™' to find 



tj"" — „ 
H x - a 



0(1) 0(e A ) 0(e 3 ) 
0(e A ) 0(e 3 ) 
0(e 3 ) 0(e) 



(128) 



This shows that further diagonalization requires rotations by angles of 0(e 3 ) or smaller, 
which we will neglect. 
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We have found that when £ = 0(e 2 ) and Sm^ > (normal hierarchy), we need an extra 
(2, 3)-rotation to diagonalize H, and the diagonalization matrix is UVWX, which we need 
to identify with 



U 



1 

C 23 S23 
-S23 C 23 



cia s 13 e- a 
1 
■Si 3 e u Ci 3 



Cl2 S12 
-Sl2 C12 

1 



(129) 



From Ref. [7], we know that when Sm^ > and a/^m^l = 0(e x ), identification of U with 
UVW yields 



013 

0i 



12 1 



6'l3 > 
'12 ~ i 

023 ~ 6*23 
5^0, 



(130) 



where ^3 = #13 + 0. (Note that we are considering the 5 = case.) Furthermore, in the 
range a/j^m^l = 0(e~ 1 ) we have 



0i 2 = ^sin(2M + --- = 0(e 3 ), 

e~i3 = ^i3 + --' = 0(e 2 ), 
a 



(131) 



which implies 





1 










1 










1 










1 













c 2 3 


S23 







1 










1 










c 2 3 


S23 







_ S23 


c 2 3 










1 










1 







- S23 


c 2 3 



(132) 



Since X is an (2, 3)-rotation matrix, multiplying UVW from the right with X will only lead 
to a shift in #23 • Therefore, 



013 ~ Q'l3 , 

9 12 » 9' 12 , 

023 ~ 023 + X , 
5^0. 



(133) 
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B4- Third Rotation, Sm^ < Case 



In this case, we diagonalize the (1, 2) submatrix of H'". Define the matrix Y as 



Y = 



Cjj 



1 



(134) 



where 



Then, 



c v = cos rj , s v = sin rj , tan 2t] = 



_ 2(as' 12 s 13 c 13 + &£c' 12 )s 



A 2 — X" 



(135) 



= Y^H'"Y 







(as'i 2 si3Ci3 + b£c' 12 )c <p s v 

{as' 12 s 13 c 13 + b^c'^c^Cr, 



An 



AS 



(136) 



where 



III 


\"r 2 


X' « 2 


1Y — 
in 


C 2 
A 2 C r? 


- S 2 
\"q 2 


2Y — 


C 2 


- S 2 



Ay_ , 



A Y+ ■> 



(137) 



with 



X"' 

A Y ± 



(A 2 + A'/) ± J(X> 2 - A'/) 2 + 4(as' 12 s 13 c 13 + btf 12 ) 2 s 2 



(138) 



Using 



as/ l 2 s 13 c 13 



K 
x' 



a 



5m 2 21 
2a 



-5* 

c 2 3 6m 31 



sin(26»i 2 ) s 13 ci 3 = aO(e A ) , 



-aO(e 2 ) , 



c\ 2 5m 2 2l 



2 = a 0(e) , 
= aO(e 3 ) , 



(139) 



and Scf, ~ 1, we find 



tan 277 = 



2(as' 12 si 3 ci 3 + 6f ci 2 )s^ _ 26£ 



A' 2 — A'/ 



5m 2 ! 



(l + 0(, 2 )} 



(140) 
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Therefore, the angle rj is given approximately by 



rj 



8ml 



(141) 



l 31 I ""'31 1 

from which we can conclude that s v = 0(e) and c v = 0(1). The eigenvalues can also be 
expanded in e and we find 



A 1Y 



A 2Y 



6 2 f 2 

5m 3i c i3 + t\- + aO(e A ) = a 0(e) , 
om 31 

6ml lC l 2 - ^ ' »^ 



+ aO(e & ) = aO(£ J ) 



(142) 



Again, the shifts are negligible. Putting everything together, we evaluate the sizes of the 



H' Y " = a 



(143) 



elements of Hv" to find 

r 0(e) 0(e 4 ) 
0(e 3 ) 0(e 3 ) 
0(e A ) 0(e 3 ) 0(1) 

which shows that further diagonalization requires rotations by angles of 0(e 3 ) or smaller, 
which we will neglect. 

Thus, we have found that when £ = 0(e 2 ) and 5m 31 < 0, the diagonalization of H requires 
an extra (1, 2) rotation, and the diagonalization matrix is UVWY . As in the Sm^ > case, 
UVWY must be identified with 



U = 



1 

C 23 S23 
-S 2 3 C 23 



c i3 u si 3 e 
1 
-^ise" c 13 



C12 S12 
-Sl2 C12 

1 



(144) 



Again, from Ref. [7], we know that the identification of UVW with U yields 



'13 



'12 



'23 



013 > 



ci 3 / , . 

sin(20i 2 ) , 



13 



2a 



a "<P 
"23 ~ 



o, 



-13 



<5m 



2a 



1 " sin(2^ 12 ) , 



(145) 



where c' 13 = cos6>' 13 . Since F is an (1, 2)-rotation matrix, multiplying UVW from the right 
by Y will only lead to a shift in 6 12 . Therefore, 



e 



13 



13 > 
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Neutrino 



Anti-neutrino 



Sm^ > (normal hierarchy) 



6*i2 is shifted by — 



2 Sm^ 



6*23 is shifted by + 



2 Sm^ 



< (inverted hierarchy) 



6>23 is shifted by + 



2 1 (5m 2 



6*i2 is shifted by — 
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TABLE I: Matter effects from neutral current universality violation. The parameter £, defined in 
Eq. (18), gives the size of the violation. 



012 

6 



23 ~ 
5 ~ 



Cl3 


( 5m 2 21 


^13 


V 2a 


#23 






c 13 v 


. 





sin(2# 12 ) + 77 , 
sin(2# 12 ) , 



l 21 



2a 



(146) 



C. Summary of Anti-Neutrino Case 

To summarize, in contrast to the neutrino case, the main effect of including the terms 
in the effective Hamiltonian for the anti-neutrinos is to shift #23 in the Sm^ > (normal 
hierarchy) case, and 9 12 in the Sm^ < (inverted hierarchy) case. The mixing angle that 
is affected depending on the sign of Sm^ is the exact opposite of the neutrino case. In the 
Sm^ > case, the shift in #23 is given by 



im^ 2 5m 

while for the Sm^ < case, the shift in 9 12 is given by 

n ~ b A _ b A _ gg (148 ) 

Listing these results together with those for the neutrino case from the previous section, 
we obtain Table I. The accuracy of our approximation will be demonstrated later by com- 
paring our conclusions with the exact numerical results. Let us now investigate how these 
shifts in the effective mixing angles affect the oscillation probabilities. 
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V. THE OSCILLATION PROBABILITIES 



The oscillation probability from neutrino flavor v a to neutrino flavor v$ in vacuum is 
given by 



2 A 21 a \tt 12 \ TT ,2\ • 2 A 31 



P (z/q _^ = i _ 4 |[/ q2 |2 (! _ |^ q2 |2) sin 2 _ 4 ^2 ^ _ gin ^ 

+ 1\U a2 \ 2 \U aZ \ 2 ^4sin 2 ^sin 2 ^ + sinA 21 sinA 31 ^ , (149) 

for the a — (3 case, and 

%^^) = 4|f/ Q2 | 2 |f/ /32 | 2 sin 2 ^ + 4|f/ a3 r|f/ / 33rsin 2 ^ i 

+2 M(U* 3 Up 3 U a2 U; 2 ) (a sin 2 sin 2 ^ + sin A 21 sin A 3 i 

+4 J( a ,/3) (^sin 2 ^y- sin A 31 - sin 2 sin A 21 ^j , 
for the a ^ f3 case, where J( a ,/3) is the Jarskog invariant, 



(150) 



J( a ,p) = +%(u* al u pi u a2 u; 2 ) = +%(u* 2 u P2 u a3 u; 3 ) = +^(u* 3 u^u al u; 1 ) 
= -sWaUfnU^uh) = -^(u* a iUpiU a3 u; 3 ) = -%u* a3 u P3 u a2 u; 2 ) 

= -J{P, a ) , (151) 



and 



5m% (5mf 4 /eV 2 ) „ „ „ 

A,, = -^fL = 2.534 V (E J GeV) ; (L /km) , 5m 2 = m 2 - m) . (152) 

The oscillation probabilities for anti-neutrinos can be obtained by replacing U with its 
complex conjugate, which amounts to flipping the sign of the CP violating phase 8. 
The oscillation probabilities in matter are obtained by making the replacements 

U a i —> U a i , Aij — > Aij = — — L , (153) 



for the neutrinos, and 



U ai - U m , ^ - Aij = , (154) 



for the anti-neutrinos. For instance, the and survival probabilities in matter are given 
by 



PW - v v) = 1 - 4 l^l 2 (l " IM 2 ) sin 2 - 4 |f/ M3 | 2 (l - |^ 3 | 2 ) 



A 2 i . ,o /_ , ,o\ . 2 ^31 

sm — — 
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+ 2|f/ M2 | 2 |^ 3 | 2 ( 4sin2 ^ isin2 ^ 1 + sinA 2iS inA 3i ) , (155) 



- ^) = 1 - 4 l^l 2 (l - l^l 2 ) sin 2 - 4 1 1?^ 3 1 2 (l - IM 2 ) sin 



2 ^31 

2 



+ 2|C/ /l2 | 2 |C/ M3 n 4sin 2 ^sin 2 ^ + sinA 21 sinA 31 ) . (156) 

Let us calculate these probabilities in the range a/IJmjj-J = 0(e~ 1 ), using the results of the 
previous sections. 

A. Neutrino Oscillations 

Al. <5m 31 > Case 

From Eq. (155), we note that we need U^ 2 and U ^ to calculate P(z^ — > i^). When 
5 = 0, these are given by 

U n2 — C\2 C 23 — s 12 s 13 s 23 i 

= C13S23 • (157) 

For the 5m 31 > case, the effective mixing angles in the region a/l&njjjj = 0(e~ r ) are well 
approximated by [7] 

^H^)^ = !-<*->. 

n 7T C13 / 0m 21 \ . , „ , 7T 



2 cl3 V2a/ m(2M+ *= 2-° (£) ' 

^23 ^ ^23 + 4^" i 6 -^ 1 ) sin(2^i 2 ) = ^23 + 0(e) . (158) 
c 13 V 2a y 

Using S13 = 1 — 0(e 4 ), Ci 3 = 0(e 2 ), we find 

^2 ~ C12C23 - S12S23 = cos(^i2 + 6*23) , 

« . (159) 



Therefore, 



Note that 



2ro/2 , 2 \ 1 „:„2 A 21 



-> « 1 - sin 2 {2(^i2 + ^23)} sin 



7T 



,2 
21 



012 + ^23 ~ ^ + ^ 23 + ( d ) { to I Sill(2 ^ 12) + X ■ (160) 
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Using s' l3 = 1 — 0(e 4 ), c' l3 = 0(e 2 ), we find 

s = sm(6' 13 - 13 ) = s' 13 c 13 + c' 13 s 13 « ci 3 , (161) 

which shows that the 0(e) terms in 6> 12 and 9 23 other than x cancel (this only happens for 
the 5 = case considered here) and we can approximate 

6 12 + 023 ~ | + #23 + X • (162) 



Therefore, 

-> ^) ~ 1 - sin 2 {2(^ 23 + X )} sin 2 ^ . (163) 
4,2. Smli < Case 

For the Sm 31 < case, the effective mixing angles in the region a/lSm^ = 0(e~ r ) are 
well approximated by [7] 

Sia « - (^l) 9.3 = 0( £ ' 2 ) , 

9i 2 «f-(^)sin(29 12 ) = |- ( £ 3), 

#23 « 6> 23 - • (164) 
Using si 3 = 0(e 2 ), c 13 = 1 - 0(e 4 ), s 12 = 1 - 0(5 6 ), c 12 = 0(e 3 ), we find 

^2 » , 

C/ M 3 ~ S23 • (165) 

Therefore, 

-P(^ M -> f„) ~ 1 - sin 2 (2# 23 ) sin" 



2/r>a \ ■ 2 ^31 



l-sm 2 {2(fl 23 -r/)}sin 2 — . (166) 



B. Anti-Neutrino Oscillations 



To calculate the survival probability, we need C/ ju2 and C/^3 as can be seen from 
Eq. (156). When 5 = 0, we have 

U^2 = Ci 2 C 23 - Sl 2 Sl 3 S23 , 
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U^3 = C13S23 • (167) 

For the Sm^ > case, the effective mixing angles in the region a/l^mU = 0(e~ r ) are well 
approximated by [7] 

£12 « ^sin(2£ 12 ) = 0(e 3 ), 

13 « ^0 13 = 0(e 2 ) , 
a 

023 ~ ^23 + X • (168) 

Therefore, si 2 = 0(e 3 ), ci 2 = 1 — 0(e e ), is 13 = 0(e 2 ), c i3 = 1 — 0(e 4 ), and we can 
approximate 

U^ 2 ~ C 23 , 

^3 ~ S23 , (169) 



which yields 

Pip? -> P M ) w 1 - sin 2 (2(9 23 ) sin 



2 A32 



l-sin 2 {2(^ 23 + x)}sin 2 ^. (170) 



52. Sm^ < Case 



For the Sm^ < case, the effective mixing angles in the region a/lSm^l — 0(e r ) are 
well approximated by [7] 



2 V a J 2 



012 « ^f^sin(2M+^ = 0(e), 
c 13 \ la J 

023 » 023 - |^ sin(2£ 12 ) = 6 23 + 0(e) . (171) 

Using S13 = 1 — 0(e 4 ), c?i 3 = 0(e 2 ), we can approximate 

U^ 2 ~ c 12 c 23 - Si 2 S 23 = cos(0i 2 + 23 ) , 

U,s ~ , (172) 



which yields 

P(^ -> P„) « 1 - sin 2 {2(£ 12 + 23 )} sin 2 4^ • (173) 



34 



Note that 

£12 + £23 « 023 - (^J) sm(2M + fj . (174) 

Using s' 13 = 1 - 0(e 4 ), c' 13 = 0(e 2 ), we find 

= sin (0' 13 - 6»i3 ) = s' 13 ci3 - c' 13 si 3 « C13 , (175) 

which allows us to approximate 

#12 + 023 « 023 + ?7- (176) 

(Again, the cancellation of the 0(e) terms other than fj occurs only for the 5 = case 
considered here.) Therefore, 

- ^) ~ 1 - sm 2 {2(^ 23 + 77)} sin 2 ^ ■ (177) 



C. Summary of Oscillation Probabilities 

To summarize what we have found, the and survival probabilities for the Sm^ > 
(normal hierarchy) case are given by 

P{y* ~> ^) ~ 1 - sin 2 {2(0 23 - C)} sin 2 ^ , 

P(p* - ^) ~ 1 - sin 2 {2(^ 23 + C)} sin 2 ^ , (178) 
while for the < (inverted hierarchy) case, they are given by 

where we have defined 



Therefore, though the effect of a non-zero £ appears in different effective mixing angles 
depending on the mass hierarchy, and whether the particle considered is the neutrino or the 
anti- neutrino (cf. Table I), the net effect on the and survival probabilities for all cases 
is to shift 9 2 3 in the oscillation amplitude. 



» l-sin 2 {2(£ 23 + C)}sin 2 ^i, 
9,) » 1 - sin 2 {2(# 23 - C)} sin 2 ^ , (179) 



9 1 Sm ?, I 



(180) 
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Unfortunately, this shift in 6> 23 may be difficult to observe. The current experimentally 
preferred value of sin 2 (2# atm ) sin 2 (2^23) is one, with the 90% lower limit given by [13, 14] 

sin 2 (2#23) > 0.92 . (181) 

Given the shape of the function sin 2 (2# 23 ) around #23 = 7r/4, sin 2 (2# 2 3) is insensitive to small 
shifts in #23 • Indeed, because of this, the angle 6> 23 itself is ill constrained, the above limit 
translating into 

9 23 = (0.2 ~ 0.3)tt . (182) 

However, our knowledge of the value of sin 2 (26*23) is to be improved considerably in 
the near future. The long baseline neutrino oscillation experiments MINOS [8], T2K [15], 
NOz/A [16], and others [17-20] will measure sin 2 (2# atm ) = 4|t7 M3 | 2 (l - \U^\ 2 ) = 4s 2 3 c 2 3 (l - 
s 23 c i3) from — > z/ M to better than 1%, while the reactor neutrino experiments Double- 
Chooz [21], KASKA [22], Braidwood [23], etc. [24] are expected to measure sin 2 (2^ rct ) = 
4 1 C/ e 3 1 2 ( 1 — I f^e3 1 2 ) = sin 2 (2#i3) from u e — » u e to an accuracy of ±0.01. These developments 
combined will determine sin 2 (2# 2 3) to better than 1%, albeit with a two-fold degeneracy. 
This degeneracy can be broken, in principle, by determining 4|L r /i3 | 2 |£/ e3 | 2 = S23 sin 2 (26 l i 3 ) 
from the CP non- violating part of the z/ M — > v e oscillation probability [16, 25-34]. Therefore, 
a unique and accurate value of sin 2 (26*23), together with whether 623 is larger or smaller than 
7r/4, may be known. Furthermore, if the 1-Megaton Hyper-Kamiokande (HyperK) detector 
is ever constructed, a JPARC— »HyperK long-baseline experiment will improve the limits 
even further [10]. 

Even then, if the central value of sin 2 (26 ) 23 ) is too close to one, then the shift due to ( 
will be invisible. Let us assume, for the sake of argument, that a 1% shift in sin 2 ^^) is 
detectable. Since 

sin 2 {2(^3 ± C)} = sin 2 (2^ 3 ) ± 2 sin(40 23 )C , (183) 
the shift due to ( would be visible if 

|2sin(4^23)CI > 0.01 . (184) 

The size of C for p = 4.6g/cm 3 , E = 17GeV, f = 0.025, and \5m 2 31 \ = 2.5 x 10~ 3 eV 2 , for 
instance, is 

C = 0.03. (185) 
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For this shift to be visible, we must have 

| sin(4^ 23 )| > | , (186) 

or 

sin 2 (2# 23 ) < 0.993 . (187) 

If we require a 2% shift, the limit will be sin 2 (26*23) < 0.97, and a 3% shift would require 
sin 2 (2^23) < 0.93. Therefore, whether the effect we are considering can be observed or not 
depends crucially on the value of sin 2 (26*23) • 



VI. NUMERICAL RESULTS 



The discussions up to this point were all based on approximate analytical calculations. To 
illustrate the accuracy of our analytical results, we presenting here the results of a numerical 
calculation of the effective mass-squared-differences, effective mixing angles, and oscillation 
probabilities. 

As inputs, we use the following: For #23, we consider the two cases 

sin 2 (20 23 ) = 1 , and sin 2 (26 , 23 ) = 0.92 with 6*23 < ^ . (188) 
The values of 9^ for these cases are 



6*23 = ^, and 6*23 = 0.204 tt . (189) 
For £, we compare the two cases 

£ = , and f = 0.025 . (190) 

£ = 0.025 corresponds to the central value of the CHARM/CHARM II constraint, Eq. (1). 
The remaining parameters are fixed to (see Ref. [7] and references therein) : 

5m 2 21 = 8.2 x 10~ 5 eV 2 , 

\5m 2 31 \ = 2.5 x 10" 3 eV 2 , 

tan 2 6*12 = 0.4 , 

sin 2 (20 13 ) = 0.16 , 

5 = 0, 
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FIG. 1: The effective mass-squared-differences, effective mixing angles, and oscillation probabilities 
for the case sin 2 (2(923) = 1. The other input parameters are given in Eq. (191). The £ = case is 
plotted with black dashed lines, while the £ = 0.025 case is plotted with gray solid lines. 
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FIG. 2: The effective mass-squared-differences, effective mixing angles, and oscillation probabilities 
for the case sin 2 (2#23) = 0.92 with #23 < tt/4. The other input parameters are given in Eq. (191). 
The £ = case is plotted with black dashed lines, while the £ = 0.025 case is plotted with gray- 
solid lines. 
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p = 4.6g/cm 3 , 
L = 9120 km. 



(191) 



The baseline length of L = 9120 km is the distance from Fermilab to Kamioka, Japan, and 
the mass density of p = 4.6g/cm 3 is the average mass density along this baseline calculated 
from the Preliminary Earth Reference Model [35]. 

Figs. 1 and 2 show the energy dependence of the effective mass-squared differences, effec- 
tive mixing angles, and oscillation probabilities of the neutrinos: Fig. 1 for the sin 2 (2^23) = 1 
case, and Fig. 2 for the sin 2 (2^23) = 0.92 (#23 < tt/4) case. In the figures, the £ = case is 
plotted in broken black lines while the £ = 0.025 case is plotted in solid gray lines. As can 
be clearly seen from the graphs in the top rows of both figures, the effective mass-squared 
differences are little affected by £ as expected. On the other hand, the graphs in the middle 
rows show that of the effective mixing angles, 9 12 is shifted in the negative direction when 
Sni^i > (normal hierarchy), while 623 is shifted in the positive direction when Sm^ < 
(inverted hierarchy), as tabulated in Table. I. However, the graphs on the bottom rows show 
that these shifts in the mixing angles are virtually invisible in the oscillation probabilities 
when sin (2# 23 ) = 1, but quite visible when sin 2 (2^ 23 ) = 0.92, again as expected. 

Numerical calculations for the anti-neutrino case also confirm the accuracy of our ana- 
lytical results, though we will not present them here. 

VII. FERMILAB -► HYPER-KAMIOKANDE 

If the value of sin 2 (2# 23 ) is not too close to one, then matter effects due to neutral 
current universality violation will lead to shifts in the oscillation probabilities, as shown 
in the bottom row of Fig. 2. Let us now ask whether such shifts are observable in long 
baseline neutrino oscillation experiments. In the following, we will assume that sin 2 (26*23) = 
0.92 (#23 < 7r/4), which was the value used in Fig. 2, and that it is accurately known. 

The effect we would like to see only appears at mass densities and energies at which 



5ml, I < a = 



(7.6324 x 10" 5 eV 2 ) x 



( 



g/cm 3 



P 




(192) 



or 




(193) 
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FIG. 3: The NUMI beam in its high energy mode. The vertical axis gives the number of expected 
charged current events at the MINOS site (732 km away from Fermilab) per kiloton of detector 
material, per year, per GeV bin without any oscillation. 

for l&njjjj = (2 ~ 3) x lCT 3 eV 2 . Since the mass density of the Earth's crust and mantle 
are 3 ~ 5g/cm 3 [35], this requires the beam energy to be larger than ~ 10 GeV. At these 
energies, the position of the first oscillation peak (dip) is determined by the condition 



which translates to 



or 



g/cm 2 



2E 



x 



L ~ 71 , 



- — I ~ 3 x TO 4 , 

km , 



L ~ 10 4 km . 



(194) 



(195) 



(196) 



Therefore, the experiment we need to consider is such that a neutrino beam of energy in 
excess of 10 GeV is aimed at a detector about 10,000 km away. 

At this point, we note that a beam with the required energies is already available at 
Fermilab. Fig. 3 is reproduced from the NUMI Technical Design Handbook [9] and shows 
the energy profile of the NUMI beam in its high energy mode. As we can see, the beam has 
considerable support in the 5 ~ 25 GeV range. The vertical axis is the expected number of 
charged current events at MINOS per kiloton of detector material, per year, per GeV bin 
without any oscillation. If a similar beam were aimed at a detector ~ 10 4 km away, which is 
more than 10 times the distance from Fermilab to MINOS, the flux will be attenuated by 
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at least 2 orders of magnitude from what is available at MINOS. Therefore, a megaton class 
detector would be required if the number of observed events is to be statistically significant. 

The planned Hyper-Kamiokande (HyperK) [10] is a megaton water-Chrenkov detector 
which would be at a distance of L = 9120 km from Fermilab. Aiming a NUMI-like beam 
from Fermilab at HyperK (the declination angle is 46 degrees) would provide the necessary 
energy, detector mass, and baseline length. So this is the setup we will consider. The average 
matter density along the baseline would be 4.6g/cm 3 , and the oscillation probability to be 
measured will be that shown in Fig. 2. 

In Fig. 4, we show the expected number of events at HyperK for 5 years of data taking. 
The dotted line indicates the expected numbers without any oscillation, and was obtained 
by rescaling the numbers from Fig. 3 to take into account the difference in baseline length, 
detector mass, and the number of years of data taking. The solid line indicates the expected 
number of events with oscillation taken into account for the normal hierarchy case with £ = 0. 
The input parameters are those listed in Eq. (191) with sin 2 (2# 23 ) = 0.92 (# 23 < tt/4). The 
£ = +0.025 and £ = —0.025 cases are shown with dashed, and dot-dashed lines, respectively. 
As one can see from the figure, the expected number if events is fairly large even at this 
distance, and even with oscillation. 

Fig. 5 shows a blowup of Fig. 4. Even at the oscillation dip at 17 GeV, the expected 
number of events is in the hundreds. Due to this significant statistics, the £ = and 
£ = ±0.025 cases are clearly distinguishable as indicated by the error bars on the £ = case. 
Therefore, this experiment can easily detect a violation in neutral current universality if it 
is as large as the CHARM/CHARM II central value. 

To see what kind of constraint this experiment could place on £, we calculate the \ 2 
between the £ = and the £ ^ cases, i.e. 



where iVj(£) is the expected number of events in the i-th GeV- wide bin, and plot the in- 
dependence of the x 2 (0 i n Fig. 6. We have restricted the bins that enter into x 2 (0 to the 8 
to 22 GeV range (14 bins), since that is the range in which the expected number of events 
fluctuates significantly with £. With 5 years of data taking, we can read off from the graph 
that the £ = and £ = ±0.005 cases are distinguishable at the 99% confidence level. This 
corresponds to a limit on universality violation at the 1% level, which will be comparable to 




8GeV<i<22 GeV 




Wfl-^(Q)] 2 

Ni(0) 



(197) 
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Fermilab2HyperK, 5 years 
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FIG. 4: The expected number of events at a Fermilab^HyperK experiment with 5 years of data 
taking. The dotted line indicates the number of events without any oscillation. The solid, dashed, 
and dot-dashed lines indicate the number of expected events with oscillation taken into account 
for the £ = 0, £ = +0.025, and £ = —0.025 cases, respectively. The mass hierarchy assumed was 
the normal hierarchy, sin 2 (2#23) = 0.92 (6*23 < tt/4), and the other input parameters were those 
listed in Eq. (191). 




FIG. 5: Blowup of Fig. 4. Error bars have been added for the £ = case. 
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FIG. 6: The \ 2 for the 14 bins from 8 to 22 GeV between the f = and £ / cases. 

the constraints from LEP/SLD [1-3] but completely model independent. For fewer years of 
data taking, the limits will be correspondingly weaker. 

We emphasize that the conclusions in this section are valid only for the sin 2 (2^23) = 
0.92 (#23 < tt/4) case. The closer sin 2 (20 2 3) is to one, the more difficult it will be to detect 
the presence or absence of £. 

VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have considered the matter effect on neutrino oscillations due to neutral 
current universality violation. It was shown that the effect of the violation appears domi- 
nantly as a shift in the effective value of #23 a t high energies, while the other effective mixing 
angles and effective mass-squared-differences are virtually unaffected. As a result, the effect 
will manifest itself as changes in the amplitudes of the oscillation probabilities, while the 
locations of the oscillation peaks and dips in distance/energy remain the same. However, 
since the amplitudes of the and P M survival probabilities are proportional to sin 2 (26 , 2 3), 
the shift in #23 would be difficult to detect if sin 2 (2#23) is too close to one. 

If the value of sin 2 (2# 2 3) is as small as 0.92, the current 90% lower bound, then a 5-year 
measurement of the survival spectrum by a Fermilab^HyperK experiment could place a 
model-independent constraint on neutral current universality violation at the 1% level. This 
would be competitive with the model-dependent constraints extracted from LEP/SLD data 
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[1-3]. 

The analysis in this paper was restricted to the 5 = case, in which the effective #23 was 
unaffected by charged-current interactions. For the 5^0 cases, one needs to account for the 
charged- current shift discussed in Ref. [7] , in addition to the neutral-current shift discussed 
in this paper, making the analysis somewhat more complicated. However, for the neutrino 
case with inverted hierarchy (Sm^ < 0), and the anti- neutrino case with normal hierarchy 
(Smli > 0), charged-current effects are always absent from #23, regardless of the value of 8. 
Therefore, using neutrinos if the hierarchy is inverted, and anti-neutrinos if the hierarchy is 
normal, can potentially provide a clean signal. 
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